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We study a slow relaxation process in a frustrated spin system in which a type of screen- 
fC I ing effect due to a frustrated environment plays an important role. This screening effect is 

D \ attributed to the highly degenerate configurations of the frustrated environment. This slow 



relaxation is due to an entropy effect and is different from those due to the energy barrier 
observed in systems such as random ferromagnets. In the present system, even if there is no 
"o^ \ energy gap, the slow relaxation still takes place. Thus, we call this phenomenon "entropic 

^ I slowing down" . Here, we study the mechanism of entropic slowing down quantitatively in 

^ I an Ising spin model with the so-called decorated bonds. The spins included in decorated 

_L I bonds (decoration spins) cause a peculiar density of states, which causes on entropy-induced 

^H \ screening effect. We analytically estimate the time scale of the system that increases ex- 

Q I ponentially with the number of decoration spins. We demonstrate the scaling of relaxation 

processes using the time scale at all temperatures including the critical point. 

^ ■ KEYWORDS: slow relaxation, Ising spin system, frustration, decorated bond system, entropy 



effect 



1. Introduction 



Q \ The mechanisms of slow relaxation in strongly interacting systems have been studied over 

the past several decades. One of them is the critical slowing down.^~^ When the correlation 

S^ ■ of an order parameter develops near a critical point, the system shows a slow relaxation. 

C^ ■ If the interactions of a system are nonuniform, as in the diluted ferromagnetic model, the 

model shows a slow relaxation even in the off-critical region. After quenching from a random 
state the order develops locally, but the domain walls between them are pinned in regions of 
relatively weak interaction. This pinning causes a very slow relaxation that may be called a 
frozen state. ^~^ Furthermore, when the system has frustration, as in spin glass, the frustration 
causes a random distribution of effective interactions and the relaxation also becomes slow.^~^^ 
In a frustrated system, many competing configurations are degenerate, and the system has 
a peculiar density of state that causes a temperature-dependent ordering structure, ^^~^^ and 
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even reentrant phase transitions.^^' ^^'^^ 

Recently, it has been found that an extremely slow relaxation takes place in a decorated 
lattice. ^^ The static properties of the decoration lattice are described by an effective coupling 
due to a decorated bond, and are the same as those of a regular lattice that has a coupling 
constant, but the dynamics of the systems are very different because the motion of decoration 
spins prevents the smooth evolution of ordering. This screening effect comes from the distri- 
bution of degrees of freedom of decoration spins, and we call it "entropic slowing down." This 
slow relaxation occurs even if there is no energy barrier. 

The purpose of the present study is to clarify the feature of the entropic slowing down, and 
to quantitatively analyze the time scale of the slowing down. Using the time scale, we demon- 
strate a scaling analysis of the relaxation processes obtained by a Monte Carlo simulation 
with the Glauber dynamics. ^^ 

Using the present slow-relaxation mechanism, we propose a new memory system that 
we call the "spin blackboard". In such a system, we can memorize all configurations at low 
temperatures and erase them by increasing the temperature. 

In § 2, we first introduce the model and present its equilibrium and dynamical properties. 
We present the probability distribution of the frustrated local structure in § 3. The effective 
time scale can be calculated analytically. In § 4, we study the relaxation of magnetization on 
a square lattice system by Monte Carlo simulation and perform the scaling plot using the 
effective time obtained in § 3. We provide a summary of this study in § 5. In Appendix, we 
analyze the flip probability by another energy unit. 

2. Model and Equilibrium Properties 

We have studied the ordering process of a decorated bond system. A decorated bond 
system consists of the so-called skeleton lattice, consisting of "system spins" and the decorated 
bonds between them. In each decorated bond, the system spins are connected by a bunch of 
paths that are mutually frustrated, as depicted in Fig. 1(a). To study the entropy effects 
on relaxation phenomena, we introduce the two-dimensional square lattice Ising spin system 
with decorated bonds depicted in Fig. 1(b). In Fig. 1(a), the frustrated local structure of 
the decorated bond unit is depicted. The circles denote "system spins" which form the square 
lattice as depicted in Fig. 1(b). The triangles denote "decoration spins" {si}, that are placed in 
the structure of a decorated bond. There are N^^ decoration spins between a nearest-neighbor 
pair of the system spins ai and cr2. We adopt the following structure of the decorated bond. 
We set half of the paths in a decorated bond in a ferromagnetic way, i.e., by two ferromagnetic 
bonds (solid lines) with a decoration spin between them. The other half of the paths are set in 
an antiferromagnetic way with ferromagnetic and antiferromagnetic bonds (dotted lines) with 
a decoration spin between them. The magnitudes of these bonds are set to be the same (J). 
As we will show below, the contributions of the ferromagnetic and antiferromagnetic paths 
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cancel out, and the effective interaction between the system spins due to the A^d paths is zero. 
In order to study the phase transition, in addition, we set the extra bond (a wavy hne) 
that causes the interaction between the system spins to be nonzero. We impose the periodic 
boundary condition on the lattice depicted in Fig. 1(b). 




(a) 



(b) 



Fig. 1. (a) Frustrated decorated bond unit. The circles and triangles denote the system and decoration 
spins, respectively. The solid, dotted, and wavy lines denote the magnetic interactions of J, —J, 
and J', respectively, (b) Two-dimensional square lattice with decorated bonds in the case of 



The Hamiltonian of a decorated bond depicted in Fig. 1(a) is 



Nd/2 TVd 

jVicr2 - J "^ Si{ai + a2) - J ^ Si (-CTi + a2) 

i=l i=Nd/2+l 



(1) 



where ai and a2 denotes the system spins and 5^ (i = 1, • • • , A/^) denote the decoration spins. 
The solid, dashed, and wavy lines in Fig. 1 denote the magnetic interactions of J, —J, and J\ 
respectively. Hereafter, we take J as the unit of energy. The effective coupling Keff between 
the system spins ai and ^2 at a temperature T — 1//3 is defined by tracing out the decoration 
spins {si}: 



where 



{s^=±l} 



A 



(4cosh2/3J)^d/2 
(3J'. 



(3) 
(4) 



Since the contributions of the left and right-half paths in Fig. 1(a) cancel out, the effective 
coupling of this system comes only from J^ The relation between the correlation function 
{(Jia2) and the effective coupling K^^ is given by 



{(Tia2) — tanhi^eff- 



(5) 
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Fig. 2. Relaxation process of the system magnetization on a square lattice system with decorated 
bonds, where the number of system spins is A^ = 20^, at T = 3, which is above the critical 
temperature. The circles indicate the relaxation process of the regular system (A^d = 0) and the 
crosses represent that of the decorated bond system of A/'d = 200. 

If J' > 0, (cricr2) is positive. If the effective coupling Keff is larger than the critical value Kc, 
the system spins have a ferromagnetic long-range order in the equilibrium state. Here, we 
consider the case of J^ = 1. 

Next, we consider the square lattice system with decorated bonds depicted in Fig. 1(b). We 
perform a single spin flip of the heat bath method of Monte Carlo simulation to consider the 
dynamical aspects in the system. We study a relaxation process of the system magnetization 

N 



^E-- (6) 



' N 

from an ordered state in which all the spins are aligned in the same direction. Here, we set the 
temperature in the paramagnetic region T = 3, i.e., Keff — 1/3 < K^— W log (l + a/2)] = 
(2.27- ••)" . At this temperature, we expect a fast relaxation to Mg = in the regular 
lattice (A/d = 0). However, in the decoration bond system {N^ > 1), we expect the effect of 
the frustrated configuration although the thermodynamic properties, such as the correlation 
functions of the system spin, are the same in both systems. 

We performed a Monte Carlo simulation of the square lattice system depicted in Fig. 1(b) 
with A^ = 20^ system spins. In Fig. 2, we compare the relaxations of Mg in the regular and 
decorated lattices. The data are obtained by taking the average of one thousand samples. 
The magnetization in the regular lattice relaxes with a relaxation time Treg — 10 Monte Carlo 
Step (MCS) we see no relaxation in the decorated bond system in this time scale. In the next 
section we will focus on the microscopic mechanism of this slow relaxation in decorated bond 
systems. 
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3. Probability Distribution of Frustrated Local Structure 

In the previous section, we showed an example of the entropic slowing down by Monte 
Carlo simulation. A large number of degenerate states cause such a slowing down. In this sec- 
tion, we analyze this slow relaxation from the viewpoint of the probability distribution of the 
decoration spins depicted in Fig. 1(a). We determine the number of states for a fixed local con- 
figuration of system spins. We depict "the parallel state" of system spins, i.e., (cri,a2) = (+, +), 
and "the antiparallel state", i.e., (cri,cr2) = (+, — ) in Figs. 3(a) and (b), respectively, where the 
black, white, and gray symbols denote the sites with +, — , and internal fields, respectively. 
We denote the number of + decoration spins in the ferromagnetic paths by tti, and that in 
the antiferromagnetic paths by n. The energy of the parallel state is given by 

E^_^_^^ (m, n) = {Na - 4m) J - f. (7) 

The energy of the parallel state does not depend on n, because the energy of each antiferro- 
magnetic path does not depend on the states of decoration spins. In the same way, the energy 
of the antiparallel state is given by 

E^f_^^ (m, n) = (-TVd + 4n) J + J^ (8) 



which is independent of m. Note that the following relations are satisfied: 



p(A^d) 



-++ (^,^) 



E^_^_^\m,n) 



,(^d) f^d 



771, n 



Ei^^Um,^ 



n 



(9) 



(10) 



Because the maximum values of m and n are N^/2, the ground state energy of the parallel 
state is given by 



4'?^(f,n)=-iVaJ-J' 



(11) 



and the number of degeneracies is 2^'^''^. Similarly, the minimum energy of the antiparallel 
state is given by 



p(Wd) 



£;;- (m,o) = -iVdJ + j'. 



(12) 




(a) 



(b) 



Fig. 3. Black, white, and gray symbols denote the sites with +, — , and internal fields, respectively. 
(a) Parallel system spin case and (b) antiparallel system spin case. 
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The lowest energy of the antiparahel state (n = 0) is higher than that of the parahel state by 

We consider the probabihty distribution of the decoration spins at a temperature T. First, 
we consider the parahel case. The probability of m up spins in the N(x/2 ferromagnetic paths 
is given by 



(2cosh2/3J)^<^/2 I m 



QTfirn) ^ —^———, ^'-^ e^-, (13) 



i^i?(n)^(^^^/MQ • a^) 



and the probability of n up spins in the antiferromagnetic paths is given by 

iVd/2 \ flY^/^ 
n 

On the other hand, if the system is in the antiparallel case, the probability of m up spins 
in the ferromagnetic paths is given by 

which is equal to i?++ (m). The probability of n up spins in the antiferromagnetic paths is 
given by 

<™H = ^^!!i_f~-/^)e-^- (16) 

+ (2cosh2/3J)^<i/2 \ n J 

which is equal to Q__^ (m). From the symmetry, the following relations hold: 



)(^d) ^_A _ n(^d) f^ 



Q^r(m) = Q'^.^>{'-^-m], (17) 



Rf^\n) = R^f^\n), (18) 

Qt'^rn) = QL^^)(m), (19) 

i?i!.^)(n) = i?L^^)(^-nV (20) 

Because R_^^ {n) and Q_^^ {m) are simple binary distributions, they are independent of tem- 
perature. The distributions become sharp when the number of decorated spins increases. On 
the other hand, Q_^^ (m) and R_^^ (n) depend on temperature. At high temperatures, they 
are maximum at A/d/4 owing to the entropy effect. At low temperatures, they are maximum 
at nearly A^d/2 and 0, respectively, because they are energetically favored states for the dis- 
tributions. 

The probability of the total distributions of the parallel and antiparallel states are given 

by 

Pil'\m,n) = Qf^\m)Rf^'\n), (21) 

Pi^j\m,n) = Q'-^^^m) R^^^^n) , (22) 
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respectively. Figure 4 shows color maps of the probability distributions P_^^ {m^n)^ 
P\_^ [m^n)^ and their overlap Pj^j^ [m^n) Pj^_ [m^n) as a function of the numbers of + 
decoration spins {m and n) for T — 100, T — 10, and T = 1 in the case of N^ — 200. 

The probability function is a binary distribution of m and n at high temperatures, e.^., 
T — 100, where the overlap is large (left columns in Fig. 4). As temperature decreases, the 
probability distributions split from the center point and the overlap decreases (center columns 
of Fig. 4). At low temperatures, e.^., T = 1, the peaks of the probability distributions move to 
the edges because decoration spins with nonfrustrated paths polarize, and the overlap becomes 
very small (right columns in Fig. 4). Therefore, the transition probability between the parallel 
and antiparallel state becomes very small at low temperatures. 

Thus, if the system is initially in the antiparallel state, there is a very small probability 
of transition to the parallel state, and the time scale of the evolution of the configuration 
becomes very long. Therefore, the system cannot reach the equilibrium state within a short 
time. 

We denote the probabilities as P++ and P+_ for the parallel and antiparallel configurations 
of the system spin obtained by tracing out the degree of freedom of decoration spins. They 
are given by 









772=0 n=0 

and they agree with the thermal equilibrium probabilities. 

To estimate the effective relaxation time of this system, we calculate the flip probability of 
the central spin in the configuration shown in Fig. 5. This spin is surrounded by two up system 
spins and two down system spins. Suppose we consider a regular spin system, z.e., N^ — 0. The 
internal fields from four neighbor system spins cancel out, and the energy difference between 
the up and down states of the central spin is zero. Thus, we call the central spin a "free spin". 
The transition probability of the center spin is 1/2 in the Glauber dynamics. ^^ However, as 
we will show below, in decorated bond systems [N^ > 0), the probability of the flip becomes 
smaller than 1/2 owing to the distribution of the surrounding decoration spins. 

The internal field on the central spin in this configuration is 

h (ni, n2, ^3, 714) = 2J {N^ - ni + 77,2 - ^3 - 714) , (25) 

where ni is the number of + spins on the ferromagnetic paths in the antiparallel state of the 
system spins. Similarly, 77,2 is that on the antiferromagnetic paths in the antiparallel state, n^ 
is that on the ferromagnetic paths in the parallel state, and 77,4 is that on the antiferromagnetic 
paths in the parallel state. 



7/16 



J. Phys. Soc. Jpn. 



Full Paper 



0.007 
0.006 
0.005 
0.004 
0.003 
0.002 
0.001 





100 
75 
50 



0.025 

0.02 

0.015 

0.01 

0.005 





25 50 75 100 



100 

75 

50 

25 





(a) 



0.007 


100 


0.006 




0.005 


75 


0.004 
0.003 


fl 50 


0.002 


25 


0.001 











(b) 



0.007 
0.006 
0.005 
0.004 
0.003 
0.002 
0.001 




100 

75 

50 

25 





(c) 



0.025 

0.02 

0.015 

0.01 

0.005 





25 50 75 100 



25 50 75 100 



25 50 75 100 



(d) 





(f) 



5e-06 


100 


4e-06 


75 


3e-06 






c 50 


2e-06 




le-06 


25 










(g) 



(h) 



(i) 



Fig. 4. (Color online) We set the number of decoration spins TVd = 200. (Top row) Probability distri- 
bution function of parallel state P^_^ (m, n) as a function of the numbers of + decoration spins m 
and n in the cases of (a) T = 100, (b) T = 10, and (c) T = 1. (Middle row) Probability distribution 
function of antiparallel state P\_J^\m^n) as a function of the numbers of + decoration spins m 
and n for (d) T = 100, (e) T = 10, and (f) T = 1. (Bottom row) Overlap probability distribution 
function of parallel state P\__^ {m^ n) and probability distribution function of antiparallel state 
Pj^_ (m, n) as a function of the numbers of + decoration spins m and n in the case of (g) T = 100, 
(h) T = 10, and (i) T = 1. Note that the ranges of values are not the same at all temperatures. 



In the thermal bath transition probability, the flip probability of the center spin Pflip is 
given by 



Vf, 



ip 



Y: QT--' (ni) R^y (-2) Q^^f ^ (ns) 4r ^ (-4) 

(711,712,723,724) 



1 + exp [-2/3/1 (ni, 712, ^3, ^4)] 



(26) 
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A^d/2 



Fig. 5. Central system spin surrounded by two up system spins and two down system spins, which 
we cah a "free spin". The black, white, and gray symbols denote +, — , and disordered spins, 
respectively. If decorated bonds are absent, the flip probability of the center system spin is exactly 
1/2. However, if decorated bonds exist, the flip probability of the center system spin is less than 
1/2 owing to the entropy effect. The number of decoration spins is N^ for each decorated bond. 



I y 

(4cosh2/?J)2^\ ^ , 



exp [4/3 J {ns - 712) 




. (27) 

1 + exp [-4/3 J (TVd - ^1 + ^2 - ^3 - ^4)] 

Figure 6 shows the dependence of the flip probability of the free spin VfHp as a function of A/^ 
at T = 5, 3, 2, 1, and 0. Here, we find that the stabilization effect is significant when A/^ is 
large at low temperatures. 



Ph 




Nh 



Fig. 6. (Color online) Flip probability of the free spin Vmp as a function of the number of decoration 
spins for several temperatures. The squares, circles, triangles, and diamonds denote the data for 
T = 5, 3, 2, and 1, respectively. The low-temperature limit of the probability is shown by the 
dotted line. 
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We define the time scale of the evolution by 

Teff =Pflip"^ (28) 

It increases rapidly with N^^. Note that the low-temperature limit of Teff is finite, because the 
slowing down is caused by the entropy effect, not by energy. Namely, 

1 /i\^^ 



This limit is shown by the dotted line in Fig. 6. 



(29) 



4. Monte Carlo Simulation 

In the previous section, we calculated the effective time of a flipping spin. We study the 
relaxation process of system magnetization on the square lattice system shown in Fig. 1(b) at 
T = 3, which is above the critical temperature in this section. We set N = 20^ as the number 
of system spins. Data are obtained by taking the average over one thousand samples in cases of 
NfX = 8, 16, 24, 32, and 40. As A/^ increases, the relaxation of system magnetization becomes 
slow. We scale these relaxation processes by the time scale Tgff (Ad) (Table I) estimated using 
eq. (28), and plot them in Fig. 7. All the relaxations of system magnetization converge in a 
scaling function, which indicates that the analysis in the previous section works well. In the 
case of Ad = 200 shown in Fig. 2, effective relaxation time is estimated to be Tgff ^ 10^^ MCS. 
Therefore, we see no relaxation within 100 MCS. 

Next, we consider the slowing down at the critical temperature and below the critical 
temperature. In the present model, the critical temperature is Tc = ^log(l + \Pi)J' because 
the effective coupling is given by K^^ — /3J^ The relaxation curves at the critical temperature 
obey a power law decay, while those below the critical temperature obey an exponential decay 
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(a) r = 2 
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(b) T = T, 
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Teff 


1.77064 X 10 


1.08986 X 10^ 


6.17373 X 10^ 


3.36029 X 10^ 


1.78598 X 10^ 



(c) T = 3 

Table I. Effective relaxation times of the free spin system Teff at T = 2, T = Tc, and T = 3 as a 
function of the number of decoration spins Ad. 
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Fig. 7. (Color online) Semilog plot of the relaxation of the system magnetization at T = 3 which 
is above the critical temperature. The squares, circles, triangles, inverted triangles and diamonds 
denote the case of A^d = 8, 16, 24, 32 and 40, respectively. Relaxation curves obey the exponential 
function: M{t) = OMe'^/^-^K 

toward spontaneous magnetization. In both cases, we find that the scahng works very well, 
as shown in Fig. 8. Therefore, we conclude that this mechanism of slowing down works at all 
temperatures. Both above and below the critical temperature, the relaxation curves can be 
fitted to the exponential function such as Ms{t) = Ae~^/^ + Mgq. In the case of T = 3, we 
obtain the following parameters: A = 0.859 and r = 7.31 at Mgq = 0. In the same way, we 
obtain the following parameters in the case of T = 2: A = 0.0609 and r = 1.52 at Mgq = 0.911. 
On the other hand, at the critical point, the relaxation curve obeys the power law such as 
Ms{t) = At-^. We obtain the following parameters: A = 0.90, b = 0.0569. 

This slowing down mechanism holds the effect of the system spins as well as that in the 
case of regular system (z.e., with no decoration spins), because the slowing down behavior is 
just a local effect. 

If we set J^ = 0, then the system is in the complete paramagnetic state, and all the spin 
configurations of system spins have the same energy. However, because of the entropic slowing 
down studied above, the system can maintain any configuration for a long time. Thus, we may 
use such a system as a new type of memory system, which we would call a spin blackboard. 
If we increase the temperature, we can erase the information of spin configuration. 

5. Conclusion 

We studied the microscopic mechanism of the increase in the relaxation time of the dec- 
orated bond system. We call this slow relaxation phenomenon "entropic slowing down". The 
origin of the entropic slowing down is the large number of degenerate configurations due to 
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Fig. 8. (Color online) Double logarithmic plot of the relaxation of the system magnetization at T = Tc 
(left panel) and semilog plot of the difference between system magnetization and the equilibrium 
value M{t) — Mgq at T = 2, which is below the critical temperature (right panel). At T = 2, the 
equilibrium value of magnetization Meq = 0.911. The symbols are the same as those in Fig. 7. 
Magnetization curves obey the power function M{t) = 0.90t~^-^^^ at the critical point (left panel) 
and the exponential function M{t) = O.OGle-^/^-^^ + 0.911 at T = 2 (right panel). 

frustration and not to the energy gap between metastable and stable states. If we set an addi- 
tional interaction such as J\ the ground state is an ordered state. However, the system cannot 
reach the ground state because of the freezing effect when it is cooled from a high-temperature 
disordered state. The system is trapped in some random configuration in a temperature region 
on the order of J, and no change occurs at low temperatures. Even if there is no additional 
interaction, this entropic slowing down also appears. We estimated effective relaxation time 
as a function of temperature and the number of decoration spins analytically. We also studied 
the relaxation process by a Monte Carlo simulation by the heat bath method. The scaling 
analysis was successfully performed at all temperatures including the critical temperature. 

Originally, the above model with decoration spins has been proposed to explain the slowing 
down in the system where reentrant phase transition occurs. ^^ We expect that this entropic 
slowing down also appears in more complicated systems such as spin glass. There, some parts 
of the random spin system are less frustrated where the spins are strongly correlated, while 
other parts are highly frustrated where the spins remain disordered. We can interpret the 
former parts by "system spins" and the latter by "decoration spins". In this course-grained 
picture, the spin order can be formed in the network of the former parts. 

Moreover, we proposed a new type of memory system, i.e. the spin blackboard in which 
any spin configuration can be stored owing to the entropy effect. If we can prepare a "micelle" 
type lattice in which a system spin is surrounded by many decoration spins, a spin blackboard 
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device can be made. We expect that some examples of entropic slowing down will be realized 
and the spin blackboard behavior will be demonstrated experimentally. 

The thermal annealing method^^'^^ is adopted in many problems: it is a very efficient 
method of obtaining a stable state. However, the thermal annealing method would not be 
efficient for obtaining the ground state in systems with entropic slowing down. A number of 
researchers have studied the quantum annealing method^^"^^ of obtaining the ground state 
of random systems as a substitute for thermal annealing. We have found that the quantum 
annealing method is successful in determining the ground state of this entropic slowing down 
system. The mechanism of this method will be reported elsewhere. 
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Appendix: Another Energy Scale: N^J 

In this paper, we set the energy unit to be J. In this appendix, we analyze the flip 
probability VfHp regarding A/^ J as the energy unit because the system spins are surrounded 
A/^ decoration spins. In Fig. A-1 we depict the flip probability of the free spin as a function 
of the number of decoration spins for the reduced temperature T/N^ — O.IJ, 0.2 J, and J. 
At low reduced temperatures, a nonmonotonic behavior appears, which can be understood 
as follows. The molecular field on a free spin from the outside decoration spins is given by 
-JA^dtanh(2/3J). When a = T/ (N^J) is smah, A^d - JA^d tanh (2/3J) 2^ JA^d is smah and 
thus molecular field increases with Ad, which causes a decrease in Pfiip. On the other hand, 
when a = T/ (N^^J) is large, — JAd tanh (2/3J) ^= — 2J/a, which does not depend on Ad- 
For a given T/N^^ the temperature T increases as Ad increases and thus P^ip approaches to 
its high temperature limit P^ip = 0.5. 
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Fig. A-1. (Color online) Flip probability of the free spin Pmp as a function of the number of decoration 
spins at several temperatures. The squares, circles, and triangles denote the data for T/N^ = 
0.1 J, 0.2 J, and J, respectively. 
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